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Abstract 

The exact formula for the neutrino oscillation probability in matter with constant density, which 
was discovered by Kimura, Takamura and Yokomakura, has been applied mostly to the standard 
case with three flavor neutrino so far. In this paper applications of their formula to more general 
cases are discussed. It is shown that this formalism can be generalized to various cases where the 
matter potential have off-diagonal components, and the two non-trivial examples are given: the 
case with magnetic moments and a magnetic field and the case with non-standard interactions. It 
is pointed out that their formalism can be applied also to the case in the long baseline limit with 
matter whose density varies adiabatically as in the case of solar neutrino. 
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I. INTRODUCTION 



Neutrino oscillations in matter (See, e.g., Ref. [l| for review.) have been discussed by 
many people in the past because the oscillation probability has non-trivial behaviors in 
matter and due to the matter effect it may exhibit non-trivial enhancement which could 
be physically important. Unfortunately, it is not easy to get an analytical formula for the 
oscillation probability in the three flavor neutrino scheme in matter, and investigation of 
its behaviors has been a difficult but important problem in the phenomenology of neutrino 
oscillations. In 2002 Kimura, Takamura and Yokomakura derived a nice compact formula 0, 
3] for the neutrino oscillation probability in matter with constant density. Basically what 
they showed is that the quantity U*jUpj, which is a factor crucial to express the oscillation 
probability analytically, can be expressed as a linear combination of U*jUpj, where U a j and 
U a j stand for the matrix element of the MNS matrix in matter and in vacuum, respectively. 

However, their formula is only applicable to the standard three flavor case. In this pa- 
per we show that their result can be generalized to various cases. We also show that their 
formalism can be applied also to the case with slowly varying matter density in the limit 
of the long neutrino path. In Sect. [Til we review briefly some aspects of the oscillation 
probabilities, including a simple derivation for the formula by Kimura, Takamura and Yoko- 
makura which was given in Ref. [H, because these are used in the following sections. Their 
formalism is generalized to the various cases where the matter potential has off-diagonal 
components, and we will discuss the case with large magnetic moments and a magnetic field 
(Sect. Illlj) and the case with non-standard interactions (Sect. HV|) . In Sect. [V]we summarize 
our conclusions. 



II. GENERALITIES ABOUT OSCILLATION PROBABILITIES 

A. The case of constant density 

It has been known [B| (See also earlier works 0, 0, HI].) that after eliminating the nega- 
tive energy states by a Tani-Foldy-Wouthusen-type transformation, the Dirac equation for 
neutrinos propagating in matter is reduced to the familiar form: 



U£U~ l + A{t) ^, (1) 



where 



£ = diag (E 1 ,E 2 ,E 3 ), 
A(t) = V2G F dmg(N e (t) - N n (t)/2,-N n (t)/2,-N n (t)/2) , 

\I/ T = (u e , z/^, v T ) is the flavor eigenstate, U is the Maki-Nakagawa-Sakata (MNS) matrix, 
Ej = Jrrij + p 2 (j = 1, 2, 3) is the energy eigenvalue of each mass eigenstate, and the matter 



effect A(t) at time (or position ) t is characterized by the density N e (t) of electrons and the 
one N n (t) of neutrons, respectively. Throughout this paper we assume for simplicity that 
the density of matter is either constant or slowly varying so that its derivative is negligible. 
The 3x3 matrix on the right hand side of Eq. ([1]) can be formally diagonalized as: 

UEU- 1 + A(t) = U(t)£{t)U-\t), (2) 
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where 



£{t) = di&g fait), E 2 (t),E 3 (t 



is a diagonal matrix with the energy eigenvalues Ej(t) in the presence of the matter effect. 

First of all, let us assume that the matter density A(t) is constant. Then all the t 
dependence disappears and Eq. ([1]) can be easily solved, resulting the flavor eigenstate at 
the distance L: 

(3) 



*(L) = C/exp (-i£L) Cr 1 ^ (0). 
Thus the oscillation probability P{y a — > z/g) is given by 



P(v n 



"(3) 



Uexp(-iSL) U- 1 



all 



j3a 
a/3 -f>a/3* 



4£Re(xfX£ 

3<k 



sin 



AEj k L 



+2 Im (xfxf*) sin (j\E jk L 

j<k 



(4) 



where we have defined 



ad 



AE jk 



Ej — E k , 

and throughout this paper the indices a, (3 = (e,/x, r) and j, k = (f,2,3) stand for those 
of the flavor and mass eigenstates, respectively. Once we know the eigenvalues Ej and the 
quantity X a/3 , the oscillation probability can be expressed analytically. 1 



B. The case of adiabatically varying density 



Secondly, let us consider the case where the density of the matter varies adiabatically as 
in the case of the solar neutrino deficit phenomena. In this case, instead of Eq. Q, we get 

-i [ L s(t)dt t/ror^fo), 

Jo 



*(L) = U(L)exp 



where U(0) and U (L) stand for the effective mixing matrices at the origin t = and at the 
end point t = L. The oscillation probability is given by 

2 



P(Vn 



U(L)exp\-i / £(t)dt\U(0) 



-i 



f3a 



^U(L) Pj U(L); k U(0r a jU(0Uexp 



-i J AE(t) jk dt 



(5) 



1 In the standard case with three flavors of neutrinos in matter, the energy eigenvalues Ej can be analytically 
obtained by the root formula for a cubic equation So the only non-trivial problem in the standard 
case is to obtain the expression for , and this was done by Kimura, Takamura and Yokomakura 0, 0] • 
In general cases, however, the analytic expression for Ej is very difficult or impossible to obtain, and we 
will discuss below only examples in which the analytic expression for Ej is known. 
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Eq. (jSJ) requires in general the quantity like U (t) pjU* (t) pk which has the same flavor index 
/3 but different mass eigenstate indices j, k, and it turns out that the analytical expression 
for U (t) pjU* (t) pk is hard to obtain. However, if the length L of the neutrino path is very 
large and if | J AE(t)jkdt\ 3> 1 is satisfied for j ^ k, as in the case of the solar neutrino 
deficit phenomena, after averaging over rapid oscillations Eq. ((SJ) is reduced to 



E xf'WX'j 



e (o), 



where we have defined 



Xf a (t) = U(t) 



la.l 



In the case of the solar neutrinos deficit process v e —>■ v e during the daylight, Xj(L) at the 

end point t = L and X" Q (0) at the origin t - 
at the center of the Sun, respectively, where 



correspond to Xf in vacuum and [X 



aal 
j J© 



a/3 



a/3 







UajUpj 

fj JI* 
u aj u /3j 







are bilinear products of the elements of the mixing matrices in vacuum and at the center of 
the Sun, respectively. Thus we obtain 



E 

j 



Xf 



Xf 







Hence we see that evaluation of the quantity I™ Q in the presence of the matter effect is 
important not only in the case of constant matter density but also in the case of adiabatically 
varying density. 



C. Another derivation of the formula by Kimura, Takamura and Yokomakura 



In this subsection a systematic derivation of their formula is given because such a deriva- 
tion will be crucial for the generalizations in the following sections. 2 The arguments are 
based on the trivial identities. From the unitarity condition of the matrix U, we have 



<W = [w- 1 ] = T.Uafi; 3 = E*f- 

3 3 

Next we take the (at,/?) component of the both hand sides in Eq. 



(6) 



USU- 1 +A] = \USU~ 1 } = y U^E.UL = Y EjX. 

\ a p V \ n R ' " J P3 Z-^ J . 



a/3 



a/3 



(7) 



2 The argument here is the same as that in Ref. 
known, it is reviewed here. 



Since this derivation does not seem to be widely 
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Furthermore, we take the (a,/?) component of the square of Eq. (j2J): 



(usu- 1 + a) 2 




U£ 2 U- 1 ' 




a(3 





a/3 



2 

j 



(8) 



Putting Eqs. (EI)-® together, we have 



1 1 1 \ f Xf 
Et E 2 E 3 \\ xf 

) \xf 



rp2 Tp2 T?2 

E x h 2 h 3 



V 



[USU- 1 + A 
(USU- 1 + Af 



Sol/3 \ 



a/3 
a/3 J 



which can be easily solved by inverting the Vandermonde matrix: 



af3 



X. 

\x 



ai3 
2 

a/3 



AE 2 iAE 3 i 
AE 2i AE 32 



(E 2 E 3 , — (E 2 + E 3 ), 1 
(E 3 Ei, -(Es + Ej, 1) 
{EtE 2 , -(E 1 + E 2 ), 1) 



[USU- 1 + A 
(USU- 1 + A) 



5af3 \ 



a/3 
af3 J 



V AE 3 \AE 32 / 
[(US U- 1 + Ay] a/ 3 (j = 1, 2) on the right hand side are given by the known quantities: 

Y^E.Xf + A5 ae 5 Pe 
j 

£ E]Xf + A^E 3 (5 ae Xf + SfrXf) + A 2 5 ae 5^ 



USU' 1 + A 
(USU' 1 + A 



a/3 



aj3 



Jf3e- 



(9) 



It can be shown that Eq. (Q coincides with the original results by Kimura, Takamura and 
Yokomakura @, S] ■ 

A remark is in order on Eq. 0. Addition of a matrix cl to Eq. (J2J) where c is a constant 
and 1 is the identity matrix, or in other words, the shift 



E; 



Ej + c (j = 1,2,3), 



(10) 



should give the same result for Xf (j = 1,2,3), since Eq. (Tl0|) only affects the overall 
phase of the oscillation amplitude and the phase has to disappear in the probability. It is 
easy to show that the shift (flOl) indeed gives the same result as Eq. (Q . The proof is given 
in Appendix |A] In practical calculations below, we will always put c = —Ei, i.e., we will 
consider the mass matrix U(S — EiVjU- 1 + A instead of the original one USU- 1 + A, since 
all the diagonal elements (S — E\l)jj = AEji = Am 2 1 /2E are expressed in terms of the 
relevant variables Amjj, and therefore calculations become simpler. To save space, however, 
we will use the matrix USU- 1 + A in most of the following discussions. 



D. The case with arbitrary number of neutrinos 

It is straightforward to generalize the discussions in sect. Ill CI to the case with arbitrary 
number of neutrinos where the matter potential is diagonal in the flavor eigenstate. The 
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scheme with number of sterile neutrinos is one of the example of these cases |4|, llCj • The 
time evolution of such a scheme with TV neutrino flavors is described by 

i d ^ = (u N £ N U N l +A N ) * N , 

where ^f N = {v ai , u a2 , • ■ ■ , v aN ) is the flavor eigenstate, 

S N = diag(E 1 ,E 2 ,---,E N ) (11) 



is the energy matrix of the mass eigenstate, 

A N = diag(A x , A 2 , ■ ■ ■ , A N ) , 

is the potential matrix for the flavor eigenstate, and Un is the N x N MNS matrix. As in 
the previous sect., by taking the a, (3 components, we get 



{UnS-nUj^ 1 + An 



j 

which leads to the simultaneous equation 



aj3 



for m = 0, • • • , N — 1, 



(I 

Ex 



1 

E 2 



1 

En 



\&l &2 • ■ ■ &N 



\ Y a/3 I 



U N 8U N l + An 



a/3 



V 



(UnSUx 1 + A N 



N-l 



(12) 



K/9 / 



Eq. ( 1121) can be solved by inverting the N x N Vandermonde matrix Vn- 

\ 

a/3 



( xf \ 

xf 


= 


ya/3 





UnSnUn + -An 



\ 



U N £ N Uf + A 



N-l 



(13) 



a/3 / 



The determinant of Vn is the Vandermonde determinant TTj<fc AEjk, and therefore V 1 
can be analytically obtained as long as we know the value of Ej. The factors [(UnZ-nUn 1 + 
AnY}^ on the right hand side of Eq. (TTBl) can be expressed as functions of the energy Ej, the 



quantity X" p in vacuum and the matter potential A y , since the matrix (Un£nU n +An) 3 is a 
sum of products of the matrices [(Un^nU^Y)^ 



T.^Xf (0<£<j) and [(A 



N) m ]e V 



(A e ) m 5 eri (0 < m < j). From Eq. (fl3|) it is clear that enhancement of the oscillation 
probability due to the matter effect occurs only when some of AEj k becomes small. 
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III. THE CASE WITH LARGE MAGNETIC MOMENTS AND A MAGNETIC 
FIELD 



So far we have assumed that the potential term is diagonal in the flavor basis. We can 
generalize the present result to the cases where we have off-diagonal potential terms. One of 
such examples is the case where there are only three active neutrinos with magnetic moments 
and the magnetic field (See, e.g., Ref. jl[ for review.). The hermitian matrix 3 



M 



with 



USU- 1 B 

B ] u*£(u*y l 



B = Bn a(3 



(14) 



is the mass matrix for neutrinos and anti-neutrinos without the matter effect where neutrinos 
have the magnetic moments \x a p in the magnetic field B. Here we assume the magnetic 
interaction of Majorana type 



(15) 



and in this case the magnetic moments fi al 3 are real and anti-symmetric in flavor indices: 

Ha/3 — —f^fSa- 

If the magnetic field is constant, then the oscillation probability can be written as 



P(v A ^ Vb) = 5ab - 4 J2 ^ (xj B X K 

J<K 

+2 ]T Im (xf B X£ B * 



AB*\ ■ 2 

1 sin 



/ AEj K L 



J<K 



sin 



V 2 



(16) 



where A, B run e, /i, r, e, p,, f, and J, K run 1, • • -, 6, respectively, and Xf B = UajU^j. Ej 
( J = 1, • • • , 6) are the eigenvalues of the 6x6 matrix A4. On the other hand, if the magnetic 
field varies very slowly and if the length L of the baseline is so long that \AEjkL\ ^> 1 is 
satisfied for J ^ K, then the oscillation probability is given by 



j=i 



(17) 



Following the same arguments as before, the quantity Xf B is given by inverting the 6x6 
Vandermonde matrix Vr. 



Xf B 



V x£ B J 



K" 1 



/ s AB \ 

[M] AB 



(Mf 



(18) 



AB 



See [5( for derivation of Eq. (| 14[) from the Dirac Eq 
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As in the previous sections, \{M) j ]ab ( J = 0, • • • , 5) on the right hand side of Eq. (fT8l) can 
be expressed in terms of the known quantities X% B and Bcd, and Eqs. (|T6|) and (|T8|) are 
useful only when we know the eigenvalues Ej. 

To demonstrate the usefulness of these formulae, let us consider the case where the 
magnetic field is large at origin but is zero at the end point and the magnetic field varies 
adiabatically. For simplicity we assume that # 13 and all the CP phases vanish. 4 In this 
case the 6x6 matrix A4 in Eq. (j!4p becomes real, and we obtain the following oscillation 
probabilities: 

3 

P {y a ^ Vp ) = P {p a ^^) = J2 (U Pj ) 2 [ReU(0) aj } 2 

3=1 

P[v a - up) = P(u a - z^) = £ (U Pj ) 2 [lmU(0) a3 } 2 , (19) 

i=i 

where U(0) the 3x3 unitary matrix which diagonalizes the 3x3 matrix U£U~ l + iB(0) at 
the origin: 

U£U~ l + iB{0) = U(0)S{0)U-\0). 

In this example the energy eigenvalues are degenerate, i.e., the 6x6 energy matrix be- 
comes diag(£, £), and the oscillation probability differs from Eq. (Tl7|) because the condition 
\AEjkL\ ^> 1(J 7^ K) is not satisfied (e.g., AEjk = not only for J = K = 1 but also for 
J — 1, K = 4). Each probability in Eqs. (j!9p itself is not expressed in terms of X" a (0), but 
we find that the following relation holds: 

3 3 

P[v a - u p ) + P{y a - up) = E (%) 2 |^(0)ai| 2 = E Xf p X^(0). (20) 

j=i i=l 

Eq. fTSU]) is a new result and without the present formalism it would be hard to derive it. 
The details of derivation of Eq. f[T9l and explicit forms of X° a (0) are given in Appendix IB1 
Eq. fTSUl) may be applicable to the case where high energy astrophysical neutrinos, which are 
produced in a relatively large magnetic field, are observed on the Earth, on the assumption 
that the fluxes of neutrinos and anti-neutrinos are almost equal. 



IV. THE CASE WITH NON-STANDARD INTERACTIONS 

Another interesting application is the oscillation probability in the presence of new physics 
in propagation [ll|, EH- In this case the mass matrix is given by 

USU- 1 + A NP (21) 



4 In the presence of the magnetic interaction ([T5|l of Majorana type, the two CP phases, which are absorbed 
by redefinition of the charged lepton fields in the standard case, cannot be absorbed and therefore become 
physical. Here, however, we will assume for simplicity that these CP phases vanish. 



s 



where 



/ 1 + e ee 






e* 








e* 





The dimensionless quantities e a p stand for possible deviation from the standard matter effect. 
Also in this case the oscillation probability is given by Eqs. (jlj) and where the standard 
potential matrix A has to be replaced by -4.Arp._The extra complication compared to the 
standard case is calculations of the eigenvalues Ej and the elements [{USU^ 1 + ANp) m ]a/3 
(m = l,2). 

Again to demonstrate the usefulness of the formalism, here we will discuss for simplicity 
the case in which the eigenvalues are the roots of a quadratic equation. It is known [13] that 
the constraints on the three parameters e ee , e PT , e TT from various experimental data are weak 
and they could be as large as 0{1). In Ref. [l4| it was found that large values (~ 0{1)) of 



the parameters e ee , e eT , e TT are consistent with all the experimental data including those of 
the atmospheric neutrino data, provided that one of the eigenvalues of the matrix f[2"Tl) at 
high energy limit, i.e., Anp, becomes zero. Simplifying even further, here we will neglect the 
parameters e efJL1 e MM , e MT which are smaller than 0(1CT 2 ) and we will consider the potential 
matrix 



A 



NP 



A 



( 1 + e ee 




e er 


Ctt 



(22) 



where A = ^/2GfN £ , the three parameters e ee , e er , e TT are constrained in such a way that 
two of the three eigenvalues become zero. We will assume that N e is constant, and we 
will take the limit Am^ — > 0. The oscillation probability Piy^ —>■ v e ) in this case can be 
analytically expressed and is given by 



P(Vu 



-4Re (AT AT *) sin 

-4Re (AT AT*) sin 
8A(AE 31 ) 2 



A_L 



4Re(xH e Xr) sin 2 



(A + — A_L)L 



A + A_(A + -A_; 



x sm 



A_L 



sm 



sm 



(A 



sin(arg(e eA1 ) + 5) 
+ - A_)L" 



(23) 



Eq. (12^1) is another new result and it would be difficult to obtain it without using the 
present formalism. The details of derivation of Eq. (|2"3"|) . explanation of the notations and 
the explicit forms of all the variables in Eq. (T2"3"|) are described in Appendix O 



V. CONCLUSIONS 

The essence of the exact formula for the neutrino oscillation probability in constant 
matter which was discovered by Kimura, Takamura and Yokomakura lies in the fact that 
the combination Xf = U a W p i* of the mixing matrix elements in matter can be expressed 
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as polynomials in the same quantity Xf = U aj U Pj * in vacuum. In this paper we have 
discussed applications of their formalism to more general cases. We have pointed out that 
their formalism can be useful for the cases in matter not only with constant density but 
also with density which varies adiabatically as in the case of the solar neutrino problem, 
after taking the limit of the long neutrino path. We have shown that their formalism can be 
generalized to the cases where the matter potential has off-diagonal components. As concrete 
non-trivial examples, we discussed the case with magnetic moments and a magnetic field, 
and the case with non-standard interactions. The application of the present formalism to the 
case with unitarity violation has been discussed elsewhere [IBj]. The formalism by Kimura, 
Takamura and Yokomakura is quite general and can be applicable to many problems in 
neutrino oscillation phenomenology. 



APPENDIX A: PROOF THAT EQ. (JTDD GIVES THE SAME ([9]) 

In this appendix we show that Eq. (TTTJT) gives the same result for Xf (j = 1,2,3). The 

U(i + ci) LJ- 1 = USU- 1 + A+cl 



value of Xf (j = 1,2,3) for 



becomes at most quadratic 5 in c, and all one has to do is to show that the coefficients of the 
terms linear and quadratic in c vanish. Let us introduce the notation 



-i 



1 1 1 

Ex + c E 2 + c E 3 + c 
(E 1 + c) 2 (E 2 + c) 2 (E 3 + c) 2 

( <w \ 

[USU- 1 + A + cl 
(USU^ + A+clf 



\ 



a/3 



( F -l ) (0) + c(y -l ) (l) +c 2 (y -l ) (2) 



where is the coefficient of the inverted Vandermonde matrix which is k-th order in c, 
and Bj is the coefficient of the vector (USU^ 1 + A + clY which is k-th order in c. Then 
the terms linear in c are given by 



+ (y-i)(o)£(i) 

/ 1 



A£oiA£ 3 i 
—„ 

V AE 3 \AE 3 2 



{E 2 + E z , -2, 0) 
-(Es + Et), +2, 0) 
-(£x+£ 2 , -2, 0) 



\ 



[USU- 1 + A 
(USU- 1 + A) 2 



5 a /3 \ 



a/3 
a/3 / 



5 Notice that all the factors AEjk are invariant under the shift (fTU|) . and the only change by this shift 
comes either from the terms EjEk or from Ej + Ek in the inverse of the Vandermonde matrix (cf. Eq. 
([9])). Hence the difference by Eq. (fT0|) is at most quadratic in c. 
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+ 



1 <+E 2 E 3 , ~{E 2 + E 3 ), +1)^ 



AE21AE31 
AE 2 \AE 32 



-E 3 E U +(E 3 + E 1 ), -1) 
+ExE 2 , -(Ex + E 2 ), +1) 





2 \UEU~ 1 + A] ap 



\ AE 3 \AE 32 

and the terms quadratic in c are given by 

(y-l)(2)£(0) + (y-l)(l)£(l) + (y-l)(0)£(2) 



AE 2X AE 3 i 
AE 2 \AE 32 

\ AE31AE32 
( 1 



0, 0) 
0, 0) 



+1, 0, 0) 



5 a B \ 



[UEU- 1 + A 
(U£U~ l + Af 



a/3 
a/3 / 



+ 



AS,W* + *' - 2 ' 0) 
' (-(S3 + S1), +2, 0) 



A.E21 A^E'32 



-{E x + E 2 , -2, 0) 



2 [LOT" 1 + .A] 




+ 



V A£ 3 iA£ 32 

- - 1 . ~ (+E 2 E 3 , -(E 2 + E 3 ), 



AE 2 iAE 32 
V AE 31 AE 32 ' 



-£ 3 £i, +(£ 3 + £i), -1) 
+£i£ 2 , —(-Bi + E 2 ), +1) 




0. 



Thus X a/3 (J = 1,2,3) is independent of c, as is claimed. 



APPENDIX B: DERIVATION OF EQ. flT9j) 

The matrix (PUD can be rewritten as 



1/ 1 «1 W U£U ^ + iB 

2 il 1 USU- 1 - iB 



1 -il \ 
-il 1 ' 



so the problem of diagonalizing the 6x6 matrix (HM is reduced to diagonalizing the 3x3 
matrices USU^ 1 ± iB. Since we are assuming that #13 and all the CP phases vanish, all 
the matrix elements U a j and B a p = —Bp a are real, USU -1 ± iB can be diagonalized by a 
unitary matrix and its complex conjugate: 

USU^ + iB = usu- 1 
USU' 1 -iB = U*S{U*)-\ 

Therefore, we can diagonalize M. by a 6 x 6 unitary matrix U as 

"-Hie)*- 1 - 
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where 



U 



V2 



i -a 
-n 1 



U 

u* 



- u 

y/2\-iU 



iU* 

u* 



We note in passing that the reason why diagonalization of the 6x6 matrix is reduced to 
that of the 3x3 matrix is because the two matrices U£U~ 1 and B are real. 

On the other hand, without a magnetic field the 6x6 unitary matrix U is given by 



U = 



U 

u* 



U 
u 



where the CP phase S has dropped out because #13 = 0. From these we can integrate the 
equation of motion and we get the fields at the end point: 



¥(L) 
* C (L) 



U(L) 



-i<!> 











* Mo) -1 U 



*(0) 
(0) 



*(0) 
# c (0) 



where 



£(t) dt, 



and we have assumed that a large magnetic field exists at the origin whereas there is no 
magnetic field at the end point. Thus the oscillation probabilities for the adiabatic transition 
are give by: 

2 



P(u a -> Up) = P(u a -> Up) = 



lim 

L^oo 



Ue-'^lf- 1 + U*e 
2 



-i<E> 



a/3 



El^| 2 [Re(^) 



P(u n 



P(u n 



lim 

L^oo 



1. - 2 



a3 



Y,\U Pj \ 2 [lm(U aj ) 
3=1 



Hence we obtain the following relation: 



To get |t/ a j| 2 , we need the explicit expression for the eigenvalues and the quantity X® a in the 
presence of a magnetic field. In the following we will subtract E\ 1 from the energy matrix E 
because it will only change the phase of the oscillation amplitude. For simplicity we will put 
#13 = 0, #23 = 7r/4, and we will consider the limit Am^ — > 0. Defining AE jk = Am J 2 J ./2£ 
and 

— p — g 
i3 a/3 = 5/x a/3 = I p -r 
g r 
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we have the eigenvalue equation 

= |A1 - U{£ -E^U- 1 -iB\ 

= A 3 - A£ 31 A 2 - (p 2 + q 2 + r 2 )A + ^(p- qf. 

The three roots of the cubic equation (1B1|) are given by 

Ai = 2Rcos(p H — , A 2 = 2Rcos((f H — tt) H — , A 3 = 2Rcos((p tc) H — . 

3 3 3 3 3 

where 



(Bl) 



R = [(AE 31 /3) 2 + (p 2 + q 2 + r 2 )/3f 2 , 

if = (1/3) cos" 1 [{(A£ 31 /3) 3 + AE 31 (p 2 + q 2 + r 2 )/Q - AE 31 (p - q) 2 /A}/R 



The quantity X ao in the presence of a magnetic field is given by 



/ 



\ A 3 



AA 2 jAA 3 i 
AA21 AA32 
V AA 31 AA 32 



(A 2 A 3 , -(A 2 + A 3 ), 1) 
(A 3 A 1; -(A 3 + Ai), 1) 
(AiA 2 , -(Ai + A 2 ), 1) 



( 1 

r 2 

\ Y°" 



(B2) 



where 



I 2 



J 3 



U(S - E^U- 1 + iB] = AE 3 iXn a 

i act ' 

' (a = e) 

^A£ 31 /2 (a = fJi,T) 

{U{£ - E^U- 1 +iB} 2 

q 2 + r 2 (a = e) 

r 2 +p 2 + (AE 31 ) 2 /2 ( a = li ) . 
p 2 + q 2 + (AE 31 ) 2 /2 (a = r) 



(AE 31 ) 2 X« a - (B 



(B3) 



(B4) 



In evaluating Y^ a , we have used the facts #13 = 0, # 23 = 7r/4, AE 2 i = 0, i3 a/ g = — Bp a , and 
that U(S — EiVjU^ 1 is a symmetric matrix. Using all these results, it is straightforward to 
obtain the explicit form for P{y a — > vp) + P(z/ Q — > up) by plugging the results of Eqs. (1B2[) . 
( 1B3I) . (jB4j) into the following (although calculations are tedious): 



P(u a -> i/ e ) + P(z/ Q 



c 2 2 xr + 

r 2 ~ s 2 ~ 1 ~ 

c 12 vna , 6 12 vaa 1 — 



P(u a _> „„) + P(p a _> ^) = + + 09 = H, r) 



where s i2 = sin#i 2 , Ci 2 = cos# 



12- 
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APPENDIX C: DERIVATION OF EQ. fl23l) 



The oscillation probability fl23l) is obtained in two steps. First we will obtain the eigenval- 
ues of the matrix (12T1) with Eq. ff22]) and then we will plug the expressions for the eigenvalues 
into Eq. §§§ with A replaced by Anp given in Eq. f[2"2l . 

Let us introduce notations for 3x3 hermitian matrices: 



A, 



-i 

1 



1 

1 



A, = 



Ac 



-i 

% 

(10 


V o o -l 



At 



/0 
0-i 
\0 z 



where A2, A5 and A7 are the standard Gell-Mann matrices whereas Ao and A9 are the notations 
which are defined only in this paper. Simple calculations show that the matrix Anp in Eq. 
f l22|) can be rewritten as 



A N p = Ae^ X9 e-^ Xs 



\ 1 4~ ^ee 4~ ^tt ^ 
^0 ~ r A9 



1 



+ if I 2 



(CI) 



where 



1 2|e er | 2 
(3 = - tan 1 — 

2 1 -)- 6 ee ^rr 



7 = 2 ar s( e ^)- 



From Eq. (ICip we see that the two potentially non-zero eigenvalues A e / and A T / of the matrix 
fl22|) are given by 



A e ' 
A r ' 



.4 



1 ~\~ £ ee ~h £ rr 



± 



l + e 



ee c tt 



-/it I 



In order for this scheme to be consistent with the atmospheric neutrino data particularly at 
high energy, which are perfectly described by vacuum oscillations, X T > has to vanish [3] . In 
this case, we have 



tan (3 



1 + e ee ' 
\ r 12 



1 + e ee 
A e / = A(l + e ee ) 



1 + 



l + e e 



A(l + eg 
cos 2 /3 



Thus we have 



^jvp = Ae i7A9 e- i/3A5 diag(A e /,0,0)e i/3A5 e- 



•17 A9 



(C2) 
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If we did not have (3 and 7, Eq. (1C2I) would be the same as the standard three flavor scheme 
in matter, which was analytically worked out in Ref. 16j in the limit of Am^ — > 0. It turns 
out that, by redefining the parametrization of the MNS matrix Eq. (1020 can be also treated 
analytically in the limit of Am|i — > as was done in Ref. [16]. The mass matrix can be 
written as 



usu 



-1 



A 



e if3x 5e - h x gusu -i e ijx 9e -i(3x 5 + diag ^ Qj ^ 



Here we introduce the following two unitary matrices: 



e ip\ 5e -ij\ 9 jj 

diag(l, 1, e iargC/ -3) u" diag(e iargC/ ^, e iargl/ ^, 1), 



U' 



where U is the 3x3 MNS matrix in the standard parametrization [17( and U" was defined 
in the second line in such a way that the elements U" x , U" 2 , U" 3 be real to be consistent with 
the standard parametrization in Ref. [17] 6 - Then we have 



U£U- l +A NP = e i7A9 e- i/3A5 diag(l,l,e iarg ^3) \u"£U"- 1 + diag (A e ,, 0, 0) 



xdiag(l, 1, e" 



-iargC/^N if3X 5 -i-yXg 



(C3) 



Before proceeding further, let us obtain the expression for the three mixing angles 9'L and 
the Dirac phase 5" in [/". Since 

/ cpe'^Uex + spe^U^ cpe'^U e2 + spe^U r2 cpe~ iry U e3 + spe^U r3 
U'=\ U, 3 

\ cpe-^Urt - s^U el cpe'^U r2 - spe^Ua cpe~^U r3 - spe l ^U e3 

where cp = cos j3, sp = sin f3, we get 

Q'[ 3 = sin" 1 \U2 = sin" 1 {cpe-^U* + sp^U r3 \ 

9'{ 2 = tan-^L/yO = tan" 1 (l^e"^ + s^U T2 \l\c^U eX + s p e*U Tl 
tajT l (U'^/U^) = tan- 1 (u^cpe^U^ - sp^U e3 \ 
-arg[/ e3 = -arg (c^e - * 7 ^ + spe^U T3 ). 



°23 

8" 



As was shown in Ref. 16], in the limit Am 21 — > 0, the matrix on the right hand side of Eq. 
]3|) can be diagonalized as follows: 



U"£U"- 1 + diag (A e /, 0, 0) - #il 



e^3^ r5 „ e ^r^ 5r -i^r 2 A2 diag 0, AE 31 ) e- l ^ X2 T 5 „e- l0 ^ X5 Tj,}e- 
e ^7 T 



AO" Mr-lje 1 ' Aa, 



+ diag (A^, 0,0) 



e <3^ diag (0j 0j AEsi ) + diag (Ae , ; 0) o) 



e<3^r>e^ A5 diag (A_, 0, A+) e'^^Tj,}^^ 



„i0" A 



where IV = diagfl, 1, e" 



A = Am 31 /2E, we have used the standard parametriza- 



tion [17j U" = e ie ^i X7 Ts"e ie '^ X5 T 5 ,}e ie "2 X2 , and the eigenvalues A ± are defined by 



A± = - (AE31 + A e /) ± -V(AS 3 i cos 2^3 - \ e ,y + (AE31 sin2^ 3 ) 2 . 



6 The element U" 2 has to be also real, but it is already satisfied because U" 2 = U r i. 



15 



Having obtained the eigenvalues, by plugging these into Eq. (J9]) with A — > Anp, Ex 
A_, E 2 -»• 0, E 3 -»• A+, we obtain X^ e : 



/Xf 

vat 



-i 



A_(A H 
1 

aIaI 



A. 



-(o, 



V A + (A + - A_ 
where Y^ e are defined by 



"(0, 



-A+, 










[ 


(A + + A_) 


, i) 










-A_, 


1} . 





( 



V fie 



-Yt + A + Kf 

A_(A + - A_) 
A + + A_)KT 



A+A_ 

Yt - a_kt 

V a + (a + -a_) y 



(t/^t/- 1 + A NP 



3-1 



fie 



and are given by 



Yt = [(AE 31 ) 2 + A(l + e ee )AE 31 )Xt + AAE 31 e* eT Xf. 

Furthermore, by introducing the notations 

£ = [(AE 31 ) 2 + A(l + e ee )AE 31 ]U ti3 \U e3 \ 
7] = AAE 31 \e eT \U^ 3 U r3 
C = AE 31 U^\U e3 \, 

we can rewrite Y 2 e = (e tS and Y 3 e = £e tS + 7]e~ 2n , where 5 is the Dirac CP phase of the 
MNS matrix U, so we have 



fie 



Xf 



AT 



A_(A + - A_ 



■it + Ve-^ +8) - A + (] 



A, A_ 



-«(2 7 +5) 



-(A + + A_)C] 



„i<5 



A + (A + -A_ 



■[ll + V e-^ +5) - A_C]. 



Notice that the phase factor e J<5 in front of each Xj 46 drops out in the oscillation probability 
P{ v n ~~ * v e) because P{y^ — > u e ) is expressed in terms of Xj e X£ e *, and the oscillation 
probability (I23I) depends only on the combination 27 + 5 = arg (e eM ) + 5. 

In the present case, the matrix U is unitary and because of this three flavor unitarity all 
the T violating terms are proportional to one factor: 

2£lm(XfXr)sin(A^L) 

3<k 

= 2 Im (Xf Xf *) [sin (AE 12 L) - sin (AE 13 L) + sin (AE 23 L 



-8 Im (Xf e Xf *) sin 



( AE 2l L\ . ( AE 31 L\ . ( AE 32 V 



sin 



sin 
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This modified Jarlskog factor Im^X^X^ 6 *) in matter can be rewritten as 

""-/ie _ 1 Tr^ivi JLe v^ e *\ — sin(27 + 5) 



A+A_(A+-A_) v 6 1 ' A+A_(A + -A_) 



A + A_(A + -A_) 
This completes derivation of Eq. (1231) . 



l^erX^ X% | sin(arg(e e ^ 
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